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Abstract. 

We present a catalog of XMM-Newton and Chandra observations of gamma-ray burst (GRB) afterglows, reduced in a common 
way using the most up-to-date calibration files and software. We focus on the continuum properties of the afterglows. We 
derive the spectral and temporal decay indices for 16 bursts. We place constraints on the burst environment and geometry. 
A comparison of the fast XMM-Newton follow-up and the late Chandra observations shows a significant dilference in those 
parameters, likely produced by a transition from jet expansion taking place between two and ten days after the burst. We do 
not observe a significant shrinking of the luminosity distribution when we correct for beaming; more burst observations are 
needed to confirm this result. We also compare our results with those obtained by BeppoSAX and SWIFT; there is no strong 
discrepancy between the afterglow fluxes observed with these satellites when we carefully take into account the different median 
observation time of each observatory. 
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1. Introduction 

Discovered at the end of the 60s liKlebesadel et al.lll973h . 
Gamma-Ray Bursts (GRBs) remained mysterious objects 
for decades. The first precise localization of long GRB 



jDezalav et al.lll996l: |] 


Souveliotouet al. II1993I) provided bv 


BeppoSAX jPiroetal. 


1 19961) led to the detection of X-ray 



11997) . and finally to the measur ement of the distance of 
these objects ( iMetzger et al. Il997l) . The afterglow emission i s 
well described by the fireba ll model Ir ees & Meszaros 1992; 
|m eszaros & Rees iil997; Panaitescu et al. 1998) and provides 
powerful diagnostics of the close environment of the burst. 
In the framework of the fireball model, the afterglow spec- 
tral and temporal evolution can be used to constrain the den- 
sity profile of the environment. Many observations point out 
a possible association between GRBs and supernovae (e.g. 
IStanek et al. ,2003: Hiorth et al. 200 3). indicating that the long 
GRB progenitor may be a massive star. In such a case the 
suiTounding medium is then expected to be the stellar wind 
arising from the star (.Chevalier & Li 1999i) . However most 
of the measurements are compatible with a constant den- 
sity profile and in some cases exclude a wind profile (e.g; 



IPanaitescu & Kumar l2002l). Onlv in a few cases does the con- 
trarv hold (e.g.lGendre et al. 120043 IPiro et al. l200,4. 

A large number of datasets available for a systematic study 
in X-ray is available. Some observations were published in sin- 
gle papers, other announced in GCN, and some not published at 
all. A uniform sample of afterglows, obtained by using the lat- 
est calibration files, is particularly required when one wants to 
use the general properties of the afterglows to constrain models. 
We have initiated a re-analysis of all X-ray afterglows observed 
so far, focusing on the burst environ ment properties. Followin g 
the release of the BeppoSAX results PasQuale et al. ,l2005l) . 
we present in this article the XMM-Newton and Chandra af- 
terglow observations carried out as of October 2004. We will 
focus on the continuum properties, and we will not discuss any 
line detections within the spectra. In a future paper (Gendre 
et al., in preparation), we will expose a methodology to detect 
lines, to assess their significance and we will discuss the fea- 
tures that we have detected. The light curves and spectra will 
be made available through a local web page'. 

This article is organized as follows. In Sec.|2]we present 
the data reduction. In Sec.|3lwe list our results. In Sec. |4]we 
discuss these results and compare them with the BeppoSAX 
catalog, before providing our conclusions. 



Send offprint requests to: B. Gendre 
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2. Data reduction and analysis 

We retrieved all public data available for GRBs that occurred 
before the 1st of October 2004, from the archives of XMM- 
Newton^ and Chandra^ observations. The complete list of re- 
trieved observations is indicated in Table ^ These observa- 
tions were calibrated using the most up-to-date software avail- 
able as of October 2004. For the XMM-Newton data, we used 
the SAS version 6.0. We focused on the data from the EPIC 
(MOS and PN) instruments ( Struder et al. 2001 ; Turner et al. 
l200lh . The calibrations were done using the tasks emchain 
and epchain. For the Chandra data, we used CIAO version 
3.1 and the calibration database CALDB version 2.28. The 
data were processed using the tasks acis^jwcess.events and 
the specific tasks used to obtain calibrated grating data'*. The 
events were filtered using all provided GTI and standard fil- 
tering criteria (GRADE==[0, 2, 3, 4, 6] for Chandra data, 
(FLAG==0&&PATTERN<=4) for XMM-Newton EPIC-PN 
data, (#XMMEA_EM&&PATTERN<=12) for XMM-Newton 
EPIC-MOS). The filtered event files were checked for flaring 
background activities. We removed any period of such activity 
using a very strict condition (a flare is defined as an increase 
by a factor of at least 5 of the mean background count rate 
during the observation, and we removed the detected flare in- 
cluding the 250 seconds that precede and follow this period). 
When less than 3000 seconds of observation remain within the 
cleaned event file, we discarded the observation. The remaining 
events were used to extract spectra and light curves. 

For non-grating data and the zero order grating data, we 
extracted spectra and light curves using circle regions. The ra- 
dius of these regions were chosen to contain at least 90% of the 
counts and take into account any possible neighboring source. 
The spectral and temporal backgrounds were extracted using 
a larger circular area free of sources at the same off-axis an- 
gle (except in the case of the zero order grating data, where 
we extracted this spectrum using an annulus centered on the 
source with inner and outer radii of about 10 and 25 pixels re- 
spectively). The XM M-Newton spectral background presents 
spatial variations (see lEhle et al. II2004I) . We have taken this 
into account by choosing a position that presents a spectral 
background similar to the source region, in order to avoid 
spurious emission or absorption features. For the other orders 
of grating data, we used the regions created by the grating 
processing tasks^. When grating data were available (either 
LETG or HETG), we obtained the spectral results by simul- 
taneously fitting the zero order spectrum and the -i-l and -1 
orders added together In such a case, we have taken into ac- 
count the cross-calibration uncertainties, which are about 20% 
(see the CIAO calibration web pages), by adding to the spec- 
tral model a multiplicative constant left free to vary within the 
range 0.8-1.2. The spectral an alysis was done using XSPEC 
version 11.3.1 llArnaudlll99^ . For the temporal analysis we 
used FTOOLS version 5.3 and a customized IDL script (see 

^ http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml 
^ http://cda.harvard.edu/chaser/mainEntry.do 
'' See http://cxc.harvard.edU/ciao3.l/threads/spectra_hetgacis/ and 
|http://cxc.harvard.edu/ciao3 . l/threads/spectra_l etgacis/| 
' see the CIAO threads web pages 



be PasQuale et al. l2005l) . We binned all data to contain at least 
20 source photons within each bin in order to use the statis- 
tic. For the spectral study, we used a model of an absorbed 
power law (with a galactic absorption and a host redshifted 
absorption). When the redshift of the host was unknown, we 
used the canonical value of z=l. For those afterglows with too 
few counts detected to produce a meaningful spectral fitting, 
we used only a galactic absorbed power law with a spectral en- 
ergy index of 1 to calibrate the count to flux conversion factor. 
For the Chandra grating data, we used only the order of the 
spectra to extract the light curve, as this was found to optimize 
the signal-to-noise ratio and thus reduce the uncertainties. 

3. Observation results 

We list in Table Hour results. We indicate undetected after- 
glows (we define as afterglow a fading bright source, or a sin- 
gle source within the error box too faint to constrain its decay) 
with a 'U' flag. For these afterglows, we indicate a flux upper 
limit assuming the source was on-axis (when there is no after- 
glow detected at other wavelengths) or at the position of the 
detected optical/radio afterglow. Bursts detected but with not 
enough counts to have a good constraint on the spectral index 
are indicated with an 'F' flag. We have derived a flux for those 
afterglows using the standard model indicated in Sec. |5] For 
the remaining bursts, we indicate if the observation was done 
by XMM-Newton ('X' flag) or by Chandra in grating mode 
('G' flag) or imaging mode (T flag). In the remainder of this 
paper, all errors are quoted at the 90% confidence level (except 
when another confidence level is indicated). Fluxes are given 
unabsorbed and within the 2.0-10.0 keV X-ray band. Readers 
interested in using our data are invited to retrieve them through 
the dedicated web page'. 

There were 31 bursts followed by XMM-Newton or 
Chandra (GRB 020321 was followed by XMM-Newton and 
Chandra). We obtained good constraints on the spectral and 
temporal indexes for 14 of these. For 3 other bursts (GRB 
991216, GRB 00 1025 A, GRB 030329), the decay index was 
not constrained but we obtained a good statistical error on the 
spectral index. For 8 bursts, the observed count rates did not al- 
low us to derive good constraints on the spectral and decay in- 
dexes. For 3 bursts (GRB 01 1 130, GRB 020531, GRB 021201), 
we did not detect any convincing afterglow. We discarded 3 
burst (GRB 010220, GRB 010222, GRB 030227) observations. 
Except for the case of GRB 011030, all results reported here 
are consistent with previous publications (whi ch are hsted in 
Table In the special case of GRB 011030. lHarrison etal.l 
(2001) reported a very different flux, not consistent with the 
count rat e observed for this sou rce. An independent analysis 
made by I d'A lessi o et al. 1 (I2OO5I1 gave results similar to ours. 
We show in Fig.^the light curves (in flux unit) of those bursts 
with a good constraint on both the spectral and decay indexes. 

3.1. Rejected observations 

GRB 010220 A strong flaring background activity occurred 
during the XMM-Newton observation, which began 14.8 hours 
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Table 1. List of observations retrieved from the archives with some basic information (configuration of the instruments, instru- 
ment that detected the prompt emission). We also report the first refereed publication on the X-ray afterglow. 



CjKd 


Instrument 


archive 

UDS-IU 


Prompt 
detection 


Fluence (band) 
erg. cm"" (keV) 


Afterglow 
Detection 


Redshift 


Ket 




cnanara AL.I0-0, rlnlLr 


jyb 




2 X lU (2U-3UU) 


X, opt, rad 


1 m 
1.U2 


1 
1 


LsKd UUU21U 


L.nanaia AL.I0-0 


0O2 


BepposAA 


0.1 X 10 (40- /OO) 


A, rad 


n /I ^ 
0.840 




LrKri UUUyzo 


Cnandia AL.I0-0 


Sys, 104/ 


Ulysses 


2.2 X 10 (25-100) 


X, opt, rad 


2.066 


2 




L-iianara AL.I0-0 


lU4o 


BeppooAA 


n vy in— S//in ^^^^ 
9.25 X 10 (40-/00) 


X, opt, rad 


1 A nn 




uKd Ul iUJU 


L-nandra AL-lo-o 


o41 1, j412 


BepposAA 


1 1 Q in— 6/'/in '7n^^ 
1.13 X 10 (40-/00) 


rad 


<3 




/^TITD Al 1 1 QA 


Cnandra AL.I5-I 


2823 


Ht i b-2 










LrKri UzUlZ / 


L.nanaia AL.lo-1 


3430, 343 / 


Ht 1 1.-2 




A, lad 






r^TJTj AOAn 1 
LrKri UzUiZl 


L.nanaia AL.lo-b 


Q A -7-7 
34 / / 


BeppoaAX 


1 r\ in— 6 //in 'vnnA 
3.0 X 10 (40- /OO) 








/^TITI A'^A/1AC 

LrKB UZU4Uj 


Chandra ACl^>-o, LblLr 


2825 


Ulysses 


1 in— 1 nn\ 

3 X 10 (25-100) 


X, opt, rad 


0.69 


15 


LrKri UZU4Z/ 


L.nandra AL,lo-o 


34y3, 34y4 


BepposAA 


^ O A 1A— 7//1A "VAAA 

< Z.y X lU (4U-/UU) 


A 




3 


LrKri UZUjjl 


L.nandra AL.lo-1 


ZoZ'^, 3o/3 


Hb 1 b-2 


1 1 r\—l(Q A A\ 
1 X lU (6-4U) 








/'~' r) D ATAO 1 T 

LrKB UZUoli 


L-nandra ACl^)-5, HblG 


4364 


Hb i b-2 


in,,! A— 5 C 1 AA\ 

3.8 X 10 (25-100) 


X, opt, rad 


1.25 


4 


/^TlTD ATlAA/f 

uKB UZiUU4 


Chandra AClo-o, HblG 


/I 1 1 /I /I AA 

4381, 44Uy 


Hb i b-2 


1 T \ / 1 A— 6 /"7 /I AA\ 

5.1 X 10 ( /-4U0) 


X, opt, rad 


2.3 


4 


LrKD UZIZUI 


L-.nanara 


a7Q/1 QTQ'^ 

^ / 0^, J 1 OJ 


TPXT 


Z X lU (,ZJ-1UUJ 








LrKri U3UZZo 


Cnandra AL.I0-0 


44Z!) 


rib i b-2 


J.I X lU (3U-4UU) 


A, opt 


1 no 


3 


GRB 030328 


Chandra ACIS-S,LETG 


4432 


HETE-2 


3.0 X 10-^(30-400) 


X, opt 


1.52 




LrKB U3UjZo 


Chandra AClo-b 


QmA "^AO 1 

ivzu, jyzi 


rib 1 b-Z 


4.0 X lU (3U-4UU) 


A, opt 







LrKB UjU/Z:5 


L-nanura al-Io-o -i 


QQ1Q QAOO 

:>V1(5, jyZZ 


rlr;, i G-Z 


U.Z X lU (3U-4UU; 


A, opt 




1 C 
lb 


LjKd UJIZZU 


L.nandra al-Io-i 




rlii i ri-Z 


1 Q \^ 1 nA\ 
V.y X lU (Zj-lUU) 


A, opt 




n 
I 


GRB 040701 


Chandra ACIS-I 


4661,4662 


HETE-2 


4.5 X 10-' (2-25) 


X 


0.2146 




f^r> T) A/1 AO 1 o 

LrKB U4UolZ 


cnandra AL.io-o 


JJD4, j3oj 


UN 1 bOKAb 


Z.J X lU (ZU-ZUU) 


A 






/"^ ri o AA 1 AT C A 

LrKB UUlUZjA 


XMM-Newton 


Al TOCTAIAI 

0128530301 


Ulysses 


T T \ / 1 A— 6 C 1 AA\ 

3.2 X 10 (25-100) 






8 


GRB 010220 


XMM-Newton 


0111520301 
0111520401 


BeppoSAX 




— 


— 


8 


/^nTD Al 1T1 1 

GKB 011211 


XMM-Newton 


AAA /I 1 A 1 A 1 

0094380101 


BeppoaAA 


T -7 in— 6 //in '7nn\ 

3. / X 10 (40- /OO) 


X, opt 


2.14 


9 


GRB 020321 


XMM-Newton 


0008820401 


BeppoSAX 


3.0 X 10-« (40-700) 








GRB 020322 


XMM-Newton 


0110980301 

Al TOCTAOAI 

0128530801 


BeppoSAX 




X, opt 




10 


GRB 0302/ / 


XMM-Newton 


0128531201 


IJN 1 bCjRAb 


Q Q 1A— 7/1A '1AA\ 

5.0 X lU {IVhlKjy)) 


V /^«f 

A, opt 




1 1 
1 1 




yVIVlIVl-i^CWlUll 


0128531501 
0128531601 




1 A lU \D\J-'-T\J\J } 


A., UpL, LaXX 


u. lUO 


1 9 


GRB 031203 


XMM-Newton 


0158160201 
0163360201 


INTEGRAL 


1.5x 10-''(20-200) 


X, opt, rad 


0.105 


13 


GRB 040106 


XMM-Newton 


0158160401 


INTEGRAL 


8.3 X 10-' (20-200) 


X, opt 




14 


GRB 040223 


XMM-Newton 


0158160601 


INTEGRAL 


4.4 X 10-' (20-200) 


X 




17 


GRB 040827 


XMM-Newton 


0164570401 


INTEGRAL 




X, opt 




16 



N.B. References are : 1 : Piro et al. "i2000), 2 : Piroet al. 12001), 3 : Amati et al. 12004), 4 : Butler et al. 12003), 5 : Klose et al. "12004)76 : 
Rau et al."(2004), 7 : Melandri et al. (2005), 8 : Watson et al. (2002a), 9 : Reeves et al. (2002), 10 : Watson et al. (2002b), 11 : 
lMereghettietal..(,2003.) . 12 : .Tiengo et al. . (■200.3.') 13 : .Watson et al. ■ (.20()4\ 14 : iGendre et al..(.2(^04a') 15 : .Mirabal et al..(.2Q0.1h . 16 : 
iDe Luca et al. N2005h . 17 : lMcGlvnn et al. I J200 j) . 18 : jButler et al. l2005l) . 



after the burst. Due to our strict flare rejection criterion, we 
discarded all the observation of this burst. 



GRB 01 0222 This observation was severely piled-up. We de- 
cided not to perform data analysis and to conservatively discard 
this observation. 



GRB 030227 The observation files gave several warnings dur- 
ing the data processing by the SAS. After a check of the event 
files, we concluded that the calibration could not be correctly 
done for this observation. We conservatively discarded this ob- 
servation. 



3.2. Particular afterglows 

GRB 00021 We used the data from BeppoSAX to constrain 
the decay index. 



GRB 001 025A Two sources were detected within the burst er- 
ror box, the brightest displaying a marginal decay (5 = 1.2+3.0, 
which implies a decrease at the 66% con fidence level) while the 
dimmer does not feature any variation. IWatson et al. I ( l2002al) 
associated the bright source with the X-ray afterglow on the 
basis of this decay trend. We report on this bright source. 
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Table 2. GRB X-ray afterglow results. For each afterglow, we give the source name, the exposure duration and the net observation 
duration, together with the results of the spectral and temporal analysis. We quote the absorption (measured at the host redshift 
or at a value of z=l when the host redshift is unknown) in excess to the galactic absorption. The X-ray absorption excess upper 
limits are given at the 90 % confidence level. We also indicate the number of detected counts for the source as a proxy of the 
detection confidence level. 



Source Flag Time Exposure Net Detected Temporal Energy Spectral Excess of Flux 

name since burst duration duration count decay spectral fit absorption 

(days) (ksec) (ksec) number index (d.o.f.) (10"' cm"-) erg s"* cm"^ 



GRB 991216 


G 


1.52 


9.8 


9.8 


1305 


(0.9) 


0.7 ± 0.2 


0.94 (56) 


<9.6 


1.92 ±0.08x10-'^ 


GRB 000210 


I 


0.84 


10.0 


10.0 


533 


1.38±0.03 


0.9 ± 0.2 


1.0(15) 


4.7 ± 1.3 


2.7 ± 0.2 X 10-'^ 


GRB 000926 


I 


2.67 


10.0 


10.0 


307 


1.9 ±0.4 


0.7 ± 0.2 


1.5 (8) 


< 3.4 


1.06 ± 0.07 X 10-" 


GRB 00 1025 A 


X 


1.88 


40.8 


33.4 


911 


(1.2) 


1.8 ±0.4 


0.90(13) 


6±3 


2.0 ± 0.2 X 10-" 


GRB 011030 


I 


10.00 


50.0 


29.6 


248 


1.9 ±0.4 


0.7 


0.1 (6) 


< 39 


5 ± 1.2 X 10-'* 


GRB 01 1130 


U 


9.77 


30.2 


30.2 















< 0.6 X 10-'5 


GRB 011211 


x 


0.50 


33.6 


33.6 


3589 


2.1 ±0.3 


1.3 ±0.1 


1.20 (67) 


<4.1 


7.3 ± 0.2 X lO-" 


GRB 020127 


F 


4.05 


10.0 


10.0 


26 














2.1 ±0.5 X 10-" 






14.53 


10.1 


10.1 


13 














0.6 ±0.2 X 10-" 


GRB 020321 


F? 


0.43 


50.0 


32.0 


227 


— 


— 


— 


— 


see text 


GRB 020322 


X 


0.62 


29.0 


28.4 


4226 


1.3 ±0.3 


1.1 ±0.1 


1.00(178) 


6.5 ± 0.8 


2.11 ±0.05 X 10-" 


GRB 020405 


G 


1.68 


51.2 


51.2 


2086 


1.6 ±0.9 


I.0±0.2 


1.24(75) 


5±2 


8.9 ±0.3x10-" 


GRB 020427 


F 


9.06 


13.9 


13.9 


55 


— 


— 


— 


— 


1.0 ±0.2 X 10-" 






17.00 


14.8 


14.8 


24 


— 


— 


— 


— 


0.4 ±0.1 X 10-" 


GRB 020531 


U 


5.14 


20.1 


20.1 


— 


— 


— 


— 


— 


< 1.2 X 10-" 


GRB 020813 


G 


0.86 


78.1 


78.1 


8041 


1.4 ±0.2 


0.83 ± 0.06 


0.93 (364) 


< 2.9 


1.63 ±0.03 xlO-" 


GRB 021004 


G 


0.85 


88.1 


88.1 


2656 


1.2 ±0.3 


1.01 ± 0.07 


1.08 (112) 


< 2.8 


4.0 ±0.2x10-" 


GRB 021201 


U 


8.73 


20.0 


20.0 


— 


— 


— 


— 


— 


< 1.2 X 10-" 


GRB 030226 


I 


1.54 


40.0 


40.0 


371 


2.7 ± 1.6 


0.9 ± 0.3 


1.1(10) 


6.6!t[;5 


3.6 ± 0.2 X 10-" 


GRB 030328 


G 


0.64 


98.1 


98.1 


1575 


1.6 ±0.3 


1.1 ± 0.2 


0.79 (39) 


<4.0 


1.70 ±0.08x10-" 


GRB 030329 


X 


37.04 


244.9 


143.5 


639 


(2) 


1.0 ±0.2 


0.98 (22) 


<0.5 


1.8 ±0.2x10-" 


GRB 030528 


F 


5.97 


26.1 


26.1 


44 










7.8 ±1.3 X 10-" 






11.80 


20.3 


20.3 












< 3.2 X 10-" 


GRB 030723 


F 


2.14 


24.8 


24.8 


76 










1.4 ± 0.7 X 10-" 






12.69 


84.1 


84.1 


70 










2.3 ± 1.7 X 10'-" 


GRB 031203 


X 


0.26 


112.8 


92.2 


6261 


0.5 ±0.1 


0.8 ±0.1 


1.04(244) 


2.6 ± 0.7 


3.95 ± 0.08 X 10-" 


GRB 031220 


F 


5.62 


40.2 


40.2 


31 










5.2 ±1.3 X 10-" 






28.47 


19.7 


19.7 












< 1.5 X 10-" 


GRB 040106 


X 


0.23 


44.6 


37.0 


15738 


1.4 ± 0.1 


0.49 ± 0.04 


1.07 (515) 


<0.5 


9.1 ±0.2 X 10-" 


GRB 040223 


X 


0.20 


45.2 


14.3 


1230 


0.7 ± 0.3 


1.6 ±0.3 


1.04(43) 


(45-88) ± 16" 


2.3 ± 0.2 X 10-" 


GRB 040701 


F 


7.77 


20.0 


20.0 


68 










1.2 ±0.2 X 10-" 






16.50 


19.0 


17.5 


39 










0.8 ±0.2 X 10-" 


GRB 040812 


F 


5.04 


10.2 


10.2 


58 










1.6 ±0.3 X 10-" 






10.16 


10.2 


10.2 


42 










1.2 ±0.2 X 10-" 


GRB 040827 


X 


0.27 


52.3 


39.3 


2899 


1.4 ±0.3 


1.3 ±0.1 


1.07(116) 


8.6 ± 1.5 


1.25 ± 0.04 X 10-" 



Note : " : A high galactic density column affect this observation. See Sec. l4.5l for details. 



GRB 011211 There was a change in the attitude of the satel- 
lite during the observation. As this change of attitude is not 
fully supported by the S AS and complicates the spectral analy- 
sis (due to the change of PN CCD-chip during the re-pointing), 
we discarded the first 5 kiloseconds of the observation (before 
the change of attitude). 

GRB 020321 This burst was detected by BeppoSA X, who 
failed to detect its afterglow din't Zand et al. ll2003b . It was 
followed-up by XMM-Newton and Chandra 10.3 and 240 
hours after the burst respectively. We detected 82 sources 
within the field of view of XMM-Newton. Twelve of these 
sour ces are located in the Wide Field Camera (WFC) error box 
jGandolfi, „2,0Q2)- The subsequent Chandra observation cov- 
ered half of the field of view of XMM-Newton. One bright 
XMM-Newton source disappeared in the Chandra field of view, 
which was claimed to be the afterglow of GRB 020321 by 
lin'tZand et 311(120031) . This source has a flux of 1.4 x 10"^ + 
0.2 X 10"'"^ erg cm"^ s"' 17 hours after the burst, and an up- 
per limit of ~ 1.1 X 10"'^ erg s"' cm"- during the Chandra 
observation. On the other hand, only about 80 % of the WFC 
error box is covered by Chandra, and in particular not all the 
bright XMM-Newton sources are in the Chandra field of view 



(as also noted bv lin't Zand et al. lEooi) . Five XMM-Newton 
sources are located outside the field of view of Chandra, but 
inside the WFC error box. There were no other observations at 
other wavelenghts to confirm this likely afterglow, and we can- 
not exclude that this variable source may be not related to GRB 
020321. 



GRB 030329 There were 3 late observations of this burst. 
The analysis was complicated by a neighboringO source that 
contaminated the spectra of the afterglow. While this contam- 
ination is low and acceptable in the first and second observa- 
tions, the last observation suff'ered strong contamination that 
prevented us from analyzing the data. Due to the late obser- 
vation time, the decay index is not constr ained using orily the 
XMM-Newton data (we refer the reader to lTiengo et al. 120031 
for the early RXTE observation). 

We report in Table |2 the spectral analysis of only the 
first observation. We did not detect any excess of absorption 
within the spectra of both observations. The spectral index is 
a = 0.85 + 0.21 ixl = 0.71, 10 d.o.f.) in the second observa- 
tion, consistent with that measured in the first observation. 
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Fig. 1. Light curves of the bursts with a good constraint on both 
the spectral and temporal indexes. We have indicated in blue 
the XMM-Newton bursts, in black the Chandra Grating bursts 
and in red the Chandra Imaging bursts (see electronic version 
for color). 



Table 3. List of GRB afterglow fluxes extrapolated or interpo- 
lated 11 and 12 hours after the burst from their spectral and 
temporal parameters (assuming there is no evolution of these 
parameters). These two times are chosen to match the works 
by De Pasquale et al. (2005) and Berger et al. (200i|l about 
Beppo-SAX and SWIFT respectively. We have separated the 
Chandra Grating (top), Chandra Imaging (center) and XMM- 
Newton (bottom) samples due to their diff'erent observation 
times, and have removed GRB 030329 due to its very late 
XMM-Newton observation. 



GRB Name 



Flux 
11 hours 
(erg cm"- s"') 



GRB 020405 
GRB 020813 
GRB 021004 
GRB 030328 



7.1 ±0 
3.93 ±0 
8.4 ±0 
2.9 + 



,3 X 10-'- 
.08 X 10-'2 
.5 X 10^" 
.2x 10-" 



Flux 
12 hours 
(erg cm"^ s"') 



6.2 ± 
3.48 ± 
7.6 ± 
2.5 ± 



.3 x 10-'2 
.07 X 10-'- 
.4x 10-'^ 
.2x 10-" 



GRB 000210 
GRB 000926 
GRB 011030 
GRB 030226 



6.7 ±0 
3.1 ±0 
1.8±0 
1.35 ±0 



.3 X 10-" 
.2 X 10-" 
.5x 10-" 
.08 X 10-" 



6.0 ± 
2.7 ± 0, 

1.6 + a 

1.10±0 



.3 X 10-" 
.2 X 10-" 
.4x 10"" 
.07 X 10-" 



GRB 001025A 
GRB 011211 
GRB 020322 
GRB 031203 
GRB 040106 
GRB 040223 
GRB 040827 



1.2 ±0 
1.72 ±0 

4.2 ±0 
4.5 ±0 
9.9 ±0 

2.3 ±0 
1.40 + 



,2 X 10-" 
,05 X 10-" 
.1 X 10-" 
.1 X lO-""* 
.2x 10-" 
,2x 10-" 
.05 x 10-" 



1.1+0 
1.46 + 
3.8 + 
4.3 + 
9.0 + 
2.2 + 
1.26 + 



.2x 10-" 
.05 x 10-" 
.1 X 10-" 
.1 X 10-" 
.2x 10-" 
.2x 10-" 
.05 X 10-" 



GRB 031 203 There were 2 observations of this burst, and we 
report in Table|2only the first observation results. The spectral 
index is a = LO + 0.3 (xl = 0.70, 18 d.o.f.) for the second ob- 
servation, consistent with that measured in the first observation. 
The decay index is compatible between the two observations. 
An expanding ring of X-ray, interpreted as the dust-scattered 
echo of llie prompt emission was also detected during the first 
observation jWatson et al. l2004 . 



GRB 031220 iMelandri et"aLl ( l2005l) have reported 2 possi- 
ble candidates which have the same flux during the first TOO. 
We report only the proposed candidate, which was not detected 
during the second TOO. 



GRB 040106 The light curve of the PN instrument displays 
a deca y with an index of L4 + 0.1. We note that Moran et al.] 
( I2OO5I) reported a very high constraint on the decay index, using 
a summed light curve. We thus summed the light curves of the 
three EPIC instruments into a single light curve. We obtained a 
decay index of 5 = L38 + 0.09 (90% confi dence level). Using 
only the 1 sigma level (as in lMoran et al. l2005i) . the error drops 
to 0.04, i.e. similar to their value. 



GRB 040701 This burst was observed twice, but several 
sources within the position error bo x show va riations. We re- 
port only the candidate proposed bv iFox et al. I (i2004) . 



4. Discussion 

4.1. Comparison between afterglows 

Twenty four afterglows were securely detected (out of 28 ob- 
servations), and one observation gave a doubtful afterglow 
detection. This is a detection efficiency of 85 % (89 % if 
the doubtful detection is indeed an afterglow). We discuss 
why some afterglows (GRB 011130, GRB 020531, and GRB 
021201) were not detected in Sec. 14. 31 and propose as an expla- 
nation a jet effect. Our sample can be separated into the XMM- 
Newton, Chandra-grating and Chandra-imaging sub-samples. 
The selection criteria to perform an observation are differ- 
ent in each of these samples. In the Chandra-grating sample, 
most of the observations were triggered because of a bright 
prompt emission detected in X-ray instruments. These bursts 
were therefore selected to be bright (see Table O in order to 
produce good quality high resolution spectra. The Chandra- 
imaging sample is mostly constituted of approved TOOs. Each 
of these TOOs have different goals (e.g. accurate localization, 
high redshift burst hunt), and one may consider this sample 
as not too biased. The XMM-Newton sample is constituted of 
bursts observed in the Director's Discretionary Time. These 
bursts (except GRB 030329) where rapidly followed because 
it was possible to observe them, without a-priori information 
on the flux level. One can consider these burst observations 
as randomly chosen within all the possible observations, and 
thus not biased. In the following, we have assumed that the 
Chandra-imaging and XMM-Newton samples do not present a 
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Table 4. Mean parameters of each sub-sample from this pa- 
per. We include the decay index (6), the spectral index {a), the 
flux at 11 hours after the burst and the observation start time. 
We also r ecall the main param eters from the BeppoSAX sam- 
ple jPe P asauale et al. Il2005h . All errors are quoted at a 68% 
confidence level. 



1.2 [ 



1 .0 ■ 





Chandra 
Imaging 


Chandra 
Grating 


XMM 


BeppoSAX 


s 


2.0 + 0.3 


1.4 ±0.1 


1.2 ±0.2 


1.3 ±0.1 e 


a 


0.8 ±0.1 


1.0 ±0.1 


1.2 ±0.2 


1.2 ±0.1 


Log(Flux) 


-11.6 ±0.3 


-11.8 ±0.3 


-12.6 ±0.2 


-12.2 ±0.1 


@ 11 hours 










Obs. time 


145.4 


95.9 


30.1 


36.5 


(ksec) 











selection bias, and we have excluded from the flux comparison 
the Chandra-grating sample, as these bursts were chosen to be 
very bright (and thus present a bias toward the bright end of the 
afterglow luminosity distribution). 

We list the fluxes, interpolated or extrapolated to a common 
time of 1 1 and 12 hours, of those bursts with a good constraint 
on the decay index in Table |3l We list the mean parameters of 
each subsample in Table|?](not taking into account the doubtful 
afterglow of GRB 020321). We recaU that the XMM-Newton 
and BeppoSAX samples are not very different in the observa- 
tion time, while the Chandra sample is observed later. In both 
Tables|3and|3this discrepancy between the XMM-Newton and 
Chandra samples appears in both the fluxes and the decay in- 
dex. The higher flux extrapolated at 1 1 hours for the Chandra 
Imaging sample can be explained as due to the steeper decay 
of these bursts : the interpolation to early times is very sensi- 
tive to a significant change of the decay properties. We stress 
that the Chandra Grating data are not strongly affected by this 
effect due to the small change of the mean decay index com- 
pared to the mean XMM-Newton decay index (see Table 
We note a marginal discrepancy not significant in the spectral 
index (the errors reported in Table |3 are quoted to the Icr level 
only) between the Chandra Imaging and the other samples. We 
also note that there is no strong difference between the XMM- 
Newton and BeppoSAX samples, which indicates that there is 
no selection bias between these two samples. Taking into ac- 
count the flux level of GRB 020321, the mean XMM-Newton 
afterglow flux decrease to IQ^'^-^*"-^ erg cm"^ s"', still of the 
same order of magnitude of the BeppoSAX one. 

Berger et al. (2005) found that SWIFT X-ray afterglows 
are fainter compared to those observed by other instru- 
ments. They explain this discrepancy by a selection effect 
on the instrument that detected the burst (BeppoSAX, IPN, 
INTEGRAL, HETE-2 versus SWIFT). They support this re- 
sult using a dataset of 14 SWIFT bursts and 49 other bursts. 
To do so, they extrapolated (or interpolated) to 12 hours after 
the burst the flux quoted in published papers or in GCNs, using 
the reported decay index (when present) or a mean value of 6 
(1.3 if the quoted afterglow flux was measured at least one hour 
after the burst, or 1 otherwise), and a simple power law decay. 
We have checked if this discrepancy also appears using the lim- 
ited sample of XMM-Newton and Chandra bursts (15 bursts). 



0.? - 




13.5 -13.0 -12.5 

Log(Flux) (erg.s"\cm 



Fig. 2. Cumulative distribution of flux of the afterglows ob- 
served by XMM-Newton (dashed line) and Swift-XRT (con- 
tinuous line) 12 hours after the burst. 



In order to take into account the bias effect toward higher fluxes 
observed in the Chandra data discussed above, we removed the 
Chandra observations and retained only the 7 XMM-Newton 
bursts. We present in Fig.|2lthe comparison of these 14 SWIFT 
bursts with the 7 XMM-Newton ones. Four Swift bursts are 
observed with a flux larger than 10"'^ erg s ' cm"^, while no 
XMM-Newton bursts are observed at this level of flux. Six of 
the Swift bursts are brighter than 5 x 10"'^ erg s"' cm"^ while 
only one XMM-Newton burst is observed at this level of bright- 
ness. It thus seems, using our very limited sample, that SWIFT 
bursts are not fainter than the XMM-Newton ones. Because the 
XMM-Newton sample is limited (7 bursts), we have checked if 
this result still holds using the more complete BeppoSAX sam- 
ple. Comparing the mean afterglow flux measured by SWIFT, 
which is lO"'^** erg cm"^ s"' at 12 hours after the burst, with 
the BeppoSAX one (lO"'^'^ erg cm"^ s"', extrapolated to 12 
hours using the main parameters listed in Table |^, we do not 
observe any strong discrepancy between the two samples. 

Moreover, we can consider the effect of the early afterglow 
evolut ion o n the comparison As r eported by iBurrows et al. I 
('2005') and IChincarini et al. I fe005'), some of the early after- 
glow light curves are not compatible with a simple power law. 
For this reason, we have discarded from the 14 SWIFT burst 
sample those for which the flux at 12 hours after the burst was 
extrapolated using an observation made at less than 3 hours 
after the burst. Four SWIFT bursts were removed. Those 4 
bursts removed from the sample have the lowest extrapolated 
fluxes. We still do not observe significant discrepancies be- 
tween SWIFT and BeppoSAX bursts : the mean SWIFT flux 
value changes to 10"'^^ erg cm"^ s"' using only the reduced 
SWIFT sample. 
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-1 
Closure relationship 

Fig. 3. The closure relationships for all bursts with constraints 
on both the spectral and temporal decay indexes. We indicate 
the closures for the three cases (Jet effects, Wind model, ISM 
model) in the three panels. Vertical lines indicate the theoretical 
expected values. Afterglows observed by Chandra are located 
at the top (in red for Chandra Imaging and in black for Chandra 
Grating data), while those observed by XMM-Newton are lo- 
cated at the bottom (in blue, see electronic version for colors). 

4.2. Time evolution of the closure relationships 

The values of the decay index {5) and the spectral index {a) are 
linked through closure relationships. These relationships de- 
pend oirdieburstenviro ( Rhoads I 
1997t ISari et al. 117995 IChevalier & iflOOol). According to 
Sari et alH lll998l) . if the burst is surrounded by a medium with 
uniform density, then : 



5-1.50-= -0.5 
(5- 1.5a = 



V > Vr 



(1) 

(2) 



In the case of a surrounding medium with a non uniform 
profile decreasing like (the wind profile, iChevalier & Li I 
ll999HChevaUer et al. llQolb . we obtain : 



5-1.50-= -0.5 
(5- 1.5a = 0.5 



V > Vr 



(3) 
(4) 
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Fig. 4. Composite light curve of an X-ray afterglow. XMM- 
Newton bursts are in blue, Chandra Grating bursts in black, 
Chandra-Imaging bursts in red. Points relative to the same af- 
terglow are connected by a solid line. The two dashed purple 
lines represent the mean decay observed from XMM-Newton 
and Chandra Imaging bursts. 



of the beam is smaller than or equ al to the (time dependent) rel- 
ativistic collimation, one obtains jRhoads 1 19971) independently 
of the density profile : 



6 - 2.Qa = V > 
5 - 2.Qa =1 v<Vr 



(5) 
(6) 



One can also estimate the geometry of the burst by a set of 
closure relationships. In case of a jet, when the opening angle 



We present the closure relationships for all bursts with good 
constraints in Fig|3 The average positions for each sub-sample 
are 5 - 2.00- = -0.93 + 0.36, -0.36 + 0.17, 0.33 + 0.32 and 
-1.14 + 0.18 for XMM-Newt on, Chandra (grating), C handra 
(imaging) and BeppoS AX (see lOe Pasauale et'al l2005l for the 
BeppoS AX case) respectively when equations |5] and |6] apply. 
When equations El and |3 apply, these positions are 6 - 
\.5a = -0.39 + 0.31, 0.10 + 0.14, 0.74 + 0.32 and -0.6 + 0.2 
respectively. 

The XMM-Newton and BeppoS AX results are very similar 
On the other hand, the mean position value changes from the 
XMM-Newton sample to the Chandra samples, with two ex- 
treme positions (XMM-Newton and Chandra imaging samples) 
and one intermediate position (Chandra grating). The mean ob- 
servation start times within these three samples are 30.1, 95.9 
and 145.4 kiloseconds for the XMM-Newton, Chandra grating 
and Chandra imaging samples. We can interpret this evolution 
either as a passing through of the cooling frequency in the X- 
ray band or as a jet signature. We discuss this in the following 
two sections, but note however that this is the first time we ob- 
serve a clear evolution with time of the global afterglow prop- 
erties from late X-ray data. 
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4.3. Burst geometry 

According to lSari et al. I(ll999l) and assuming a burst located at 
z-l, the beaming angle is : 



61 = 0.166 



no 



[El 



52 



1/8 



rad 

h 



(7) 



In Eqn. «o is the density in particle cm"^, £, 52 is the 
isotropic energy in units of 10^^ erg, and ti, is the jet break 
date expressed in days after the burst. The top panel of Fig. 
13 displays the result for the jet signature, and indicates that 
we can rule out a jet signature within XMM-Newton bursts 
and cannot rule out this signature within Chandra bursts. This 
does not exclude the possibility of a coUimated fireball in the 
XMM-Newton burst cases, but simply puts a lower limit on 
the jet opening angle. Other facts indicate a possible break in 
the Chandra burst light curves. As can be noted in Table |3 the 
Chandra bursts appear brighter than the XMM-Newton ones 
when extrapolated back to half a day after the burst. Within the 
Chandra sample, there is also a trend : the brightest afterglows 
(12 hours after the bursts) are also the ones observed latest. All 
of this may indicate that the crude extrapolation using a single 
power law is not correct, and that there is a break in the light 
curve before the observation. Moreover, the position evolution 
with time of the calculated jet closure relationships (Eqn.|5]and 
|6j as previously discussed can be interpreted as a convergence 
toward the expected value in the case of a jet signature. This is 
also clear if one produces a composite light curve from all the 
bursts observed (see Fig.0}. 

To produce Fig. |3 we first rescaled the flux levels of the 
XMM-Newton bursts observed before 10^ seconds to match 
the mean XMM-Newton flux level at 11 hours (which gives 
the value 1 in Fig.|^. We then scaled the Chandra grating flux 
levels to the mean Chandra grating flux level, and used a mul- 
tiplicative factor to smoothly connect the XMM-Newton and 
Chandra grating light curves. We then rescaled the remain- 
ing bursts (one late XMM-Newton and 3 Chandra Imaging 
observations) to connect the composite light curve. One can 
clearly see an evolution from the XMM-Newton to the Chandra 
Imaging observations. The decay changes from 1.2 + 0.2 to 
2.0 + 0.3 (A5 = 0.8 + 0.5, 90% confidence level), consistent 
with a jet break (the theoretical value is A5 = 0.83 - 1.33, de- 
pending on the position of the cooling frequency and the burst 
environment) but not with a cooling break (A(J - 0.25). We 
note however that some faint events feature a shallower appar- 
ent decay (e.g. GRB 030723), not compatible with a jet break. 
While we cannot exclude the hypothesis that some of these ob- 
ject with a shallow decay are not related to GRBs (other galac- 
tic X-ray objects are variable, see e.g. Gendre et al. 2003) or 
indeed present a large beaming angle, this shallow decay may 
also indicate ^anng afterglows like GRB 970528 (Piro et al. 
1998.) or GRB 050904 ( Watson etal. 2005; Cusumanp_e_t al. 



20051) . The Ught curve resolution does not allow us to discrim- 



inate between these 3 hypothese, and in the remaining we use 
only bursts with a good light curve resolution. Doing so, we 
thus explain the evolution of the closure relationship values 
and the steep Chandra Imaging burst decays by a jet effect 
at late times. Assuming there is no selection bias within the 



XMM-Newton and Chandra Imaging samples, and taking into 
account the mean observation time (1 day after the burst and 2 
to 10 days after the burst for XMM-Newton and Chandra re- 
spectively), we can set lower and an upper limits : 

1/8 / xl/8 



0.166 



»o 

Ei,52 



<e< 0.39 



»0 



(8) 



As suming typical va lues for «o and £',,52 (0.1, 10 respec- 
tively, lFrailetal.lf200ll) . we obtain from Eqn. ^ : 0.09 rad 
< < 0.22 rad (5.1 deg < 6 < 12.6deg)j which is i n agree- 
ment with previous works (e.g. lFrail et al. I I2OOI: Ber ger et al. I 
l2003h . 

The steep decay observed late can also explain the non- 
detections. The observations that led to non detections are per- 
formed very late (see Table|5J. At that time, one should observe 
a steep decay of the afterglow due to the jet effect. This could 
prevent a detection due to the afterglow faintness at that time. 

We have looked for a standard energy emission in the X-ray 
afterglow as in Berger et al. .(.2003.) . As the XMM-Newton and 
Chandra catalog does not present enough bursts with known 
beaming angle to have meaningful results, we h ave extended 
our sample using the BeppoSAX catalog cDe Pasauale et al.1 
2005). We have computed the X-ray luminosi t y usin g the con- 
ve rsion formula quoted in iLamb & ReichartI ilOOd) and used 
bvlBer ger et"anil2003h . We have included in our sample bursts 
neither with unknown redshift, nor bursts with poorly con- 
strained spectral and temporal decay indexes. Eighteen bursts 
were included, out of which 14 have a known beaming an- 
gle. We have computed the luminosity at a common epoch of 
1 1 hours after the burst using for each burst its spectral and 
temporal decay indices (note that Berger et al. . 2003. used a 
mean spectral index and mean values for the temporal decay 
index and redshift for unknown values). The result is shown 
in Fig. |5] In the top panel, we have plotted the isotropic lu- 
minosity distribution for all the 18 bursts indicated above. In 
the middle panel, we plotted the isotropic luminosity distri- 
bution for only those bursts with a known beaming angle (14 
bursts), and in the bottom panel, we plotted the beaming cor- 
rected luminosity distribution of these 14 bursts. We have fitted 
a Gaussian distribution to each luminosity distribution. We ob- 
tained cr - 0.5 + 0.3 for the upper panel, cr = 0.5 + 0.2 for 
the middle panel a nd cr - 0.4 + 0.2 for the lower panel. The 
values obtained by iBerger etal.l l2003) were cr = 0.7 ± 0.2 
for the isotropic luminosity distribution and cr = 0.3 + 0.1 
for the beaming-corrected luminosity distribution. We observe 
the same (within error bars) distributions when one corrects 
for beaming. However, our isotrop ic luminosity distrib ution is 
narrower than the one observed bv lBerger et al. I ('2003), while 
compatible within the error bars. In fact, its width is consis- 
tent with that of the beaming-corrected distribution. So, we 
do not observe a significant shrinking of the luminosity dis- 
tribution when t aking into account th e beaming. We thus can- 
not confirm the iBer ger et al. I ll2003l) claim for a standard re- 
lease of energy in the X-ray afterglow. The addition of sev- 
eral SWIFT bursts may reduce the error bars and may help to 
reach a conclusion on that elTect. Also, we can observe in the 
top panel a bimodal distribution compatible with the results of 
lGendre&Boer.(.20()l . 
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Fig. 5. Distributions of the isotropic luminosity for all 
BeppoSAX, XMM-Newton or Chandra bursts with good con- 
straints on the spectral and temporal decay index and known 
redshift (top panel), for a sub-sample of bursts with known 
beaming angle (middle panel), and distribution of the beaming 
corrected l uminosity of t h is sub -panel (bottom panel) using the 
method of Ifiereer et al. I (l2003l) . All luminosities were calcu- 
lated at 11 hours after the burst (in the burst rest frame), using 
the observed values of spectral index, temporal decay, redshift 
and beaming angle. The low luminosity outlier in the top panel 
is GRB 031203, not plotted in the two other panels. 



4.4. Surrounding medium 

The central and bottom panels of Fig.|3display the results for 
the wind and constant density environments respectively. 

Most of the bursts can fit both of the medium classes. This 
is mainly due to the degenerations observed for 2 closure rela- 
tionships, requiring data taken below to constrain the burst 
environment. However, we can distinguish the medium type for 
one burst. 

Using the X-ray afterglow data only, we can constrain the 
burst env ironment of GRB 040106 to be a wind environment 
jGendre et al. 2004b) . Moran et al. (20Q5) ruled out this con- 
clusion on the basis of their re-analysis, stating that the closure 
was at more than 3cr of the theoretical position. They based 
their conclusions on the Icr error bars they quoted. Using our 
analysis and the 90% confidence level as derived from a study 
of the space for all the errors (which is larger than twice the 
Icr error), we confirm that the theoretical position of the clo- 
sure is within the 90% confidence level region of the observed 
closure, and thus we claim that this burst is indeed surrounded 
by a wind profile due to the burst progenitor. 



4.5. Absorption 

From our analysis, 8 bursts (of 17) present an excess of ab- 
sorption that ranges from 0.5 x 10^' to 88 x 10^' cm"^, with 
a median value of ~ 5.8 x 10^' cm"^. We present the distribu- 
tion of absorption versus the redshift in Fig. |6l GRB 040223 
is located near- the galactic plane (1=341.6138°, b=3.1940°). 
The galactic column densit y in that direction is 6 x 10^' cm"^ 
dPickev & Lockrnan Il990 l). and varies (within one degree) be- 
tween 4.46x10^' and9. 26x10^' cm"- (as one may expect, this 
increase is due to a ch ange in the distance to the galactic disk). 
ISchlegel et all i ll 998 ^ indicated that their work cannot be used 
with 1^1 < 5°; as a diagnostic, the Nh value expected from their 
map of dust is ~ 9 x 10^'. Assuming that all the absorption 
is galactic, we obtain Nh - (1.8 ± 0.3) x 10^^, which is not 
compatible with eithe r the jPickev & Lockman ( 1990) or the 
ISchlegel et all dl 9981) values. Assuming the largest observed 
absorption value within one degree, one obtains an extragalac- 
tic absorption of Nh,z=i - (4-5 + 1 .5) x 10^^, which is larger than 
the values observed in other bursts. In the remaining, we will 
conservatively use this value. If this burst is at short distances 
(z ~ 0.1), then the extragalactic absorption value becomes of 
the same order as the one observed in other bursts. 

We have investigated whether some non detections of the 
optical afterglow may be due to the absorpt ion. Assuming 
a gala ctic gas to dust law, using the work of ISchlegel et al 1 
(Il998h . the observed excess of absorption implies an extinc- 
tion in the R band between 0.27 and 25.06 magnitudes (with a 
mean of 3.23 magnitudes, excluding the value of GRB 040223 
from the calculation due to its uncertainty). This high opti- 
cal e xtinction is not agreement with optical data and lead 
iGalama fc Wiiei^ dioOII) to propose that GRBs can des troy 
dust around them. On the other hand, 'Stratta et al. I ll2004l) in- 
dicated that the galactic gas to dust law is not suited for GRB 
studies, and indicate d tha t the correct law s may be those of 
ICalzetti et al. 1 (11994 and ICalzetti 1 (ll997l) . Using these laws, 
the R extinction is still ranging from 0.13 to 1 1.08 magnitudes 
(with a mean of 1.43 magnitudes). A burst like GRB 040223, 
with an R extinction of at least 1 1 .08 magnitudes, would not be 
detectable in the optical, even at very early times. 

Bursts with no optical detection but with a posi- 
tive detection in X-rays or radio are called Dark Bursts, 
^e PasQuale et al. ( 2003) indicated that 75% of dark bursts are 
faint bursts with a failed detection of the optical afterglow due 
to its faintness. The remaining ones are those that do not fit 
the simple fireball afterglow spectra l model (unabsorbed bro- 
ken power laws, see lSari et al. 1 19981) because of a depletion in 
the optical. According to'De Pasau ale et al. ] 12003'), these truly 
dark bursts account for the remaining 25 % of dark bursts. They 
can be distant bursts (z > 5, so that the Lyo. break is redshifted 
into the optical range) or absorbed bursts (the effect of absorp- 
tion would be more severe in optical than in X-rays or radio, 
an d can prevent an optical d etection). This result was confirmed 
bv Llakobsson et ain(l2004 i) using a different method. 

Eight bursts of our sample do not present an optical af- 
terglow. GRB 040223 is located near the galactic plane, and 
the galactic absorption can explain the non-detection, as indi- 
cated before. Two of the remaining 7 dark bursts are absorbed 
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Fig. 6. X-ray absorption versus redshift. We indicate in this fig- 
ure all detected absorptions or the corresponding upper limits 
for all bursts bright enough for a spectral study. The redshift 
was assumed to be one when not available. 




Log X-ray flux (erg/s/cm/cm) 

Fig. 7. Optical versus X-ray fluxes of GRB afterglows 11 
hours after the burst. Squares and open circles repre- 
sent XMM-Newton and Beppo-SAX data (extracted from 
IPe Pasquale et al.1l2003h respectively. Lines indicate the best 
fit relationships (see text for details). 



(one of these tvvo burs ts has been classified as truly dark by 
Ijakobss on et al. 1 12004'). This represents ~ 28 % of the dark 
burst sample. We used a very limited sample, and that this re- 
sult should be considered as a trend. This number is of the same 
order as the 25 % of truly dark bursts (the counting statistic un- 
certainty is ~ 14 %). The truly dark bursts may thus be only 
absorbed bursts, whereas we cannot exclude that a small frac- 
tion (<~ 14%) of dark bursts may in fact be distant (z> 5) 
bursts. 



We present the optical versus X-ray fluxes diagram in Fig. 
Because of the bias discussed earlier and the possible jet 
effects in the light curves, we used only the XMM-Newton 
sample, together with the BeppoSAX sample, to build this fig- 
ure. All the values have been corrected for the galactic ab- 
sorption (X-ray) or e xtinction (R band), using the work of 
ISchlegel etal.l ( Il998l) . Not all XMM-Newton bursts are dis- 
played in this figure : three GRBs in our sample do not have de- 
tected optical afterglows (or the detections are very late). We do 
not include the corresponding upper limits in Fig.EJdue to the 
poor constraints they give (they are R< 18 at 11 hours after the 
burst once corrected for the galactic reddening). We have fitted 
these data points using an empirical power law. The best fitted 
relationships are indicated in the figure. To derive these rela- 
tionships, we have included (dashed line) and excluded (solid 
line) from the fit the upper limits. We can use these laws to test 
the hypothesis of absorption around faint optical afterglows. 
To do so, we focus on GRB 020322. This burst is located in the 
"dark" side of the figure : it is a normally bright X-ray after- 
glow, while it displays a faint optical afterglow. It also displays 
an excess of X-ray absorp tion (see Table [I] and Watson et^l^ 
,2002b). Using the laws of ICalzetti et al. (1994.) and Calzetti 
(1.1997), we calculated an extinction of 1.6 + 0.4 in the R band. 
The error on this value is calculated taking into account all the 
uncertainties of the spectral fit (i.e. on the galactic absorption, 
on the intrinsic absorption and on the spectral index). The best 
fitted relationships imply an extinction of 0.5 -2.5 for this burst, 
compatible with our finding from the X-ray. As other bursts 
present an extragalactic absorption higher than the one of GRB 
020322, some of the dark bursts located in the same area of this 
diagram may be absorbed bursts. Note that this does not imply 
that all dark bursts are absorbed ones. In fact, ~ 50% of bursts 
present an excess of absorption, and 5 of these are not dark : the 
presence of absorption does not conflict in all the cases with the 
optical afterglow detection. 

5. Conclusion 

We have presented a catalog of X-ray afterglows observed with 
XMM-Newton and Chandra. We have derived the absorption 
in the host galaxy, the decay and spectral indexes, and the ob- 
served unabsorbed fluxes. 

We have observed variations of the global properties of the 
X-ray afterglows within the samples that we explain by a jet 
effect in the sample. This is observed from both a variation of 
the closure relationships and the steepening of the mean de- 
cay index for bursts observed late. We have shown that jet ef- 
fects are not present within one day after the burst. Taking ad- 
vantage of the late observation time of Chandra, we find that 
the jet signature occurs between two and ten days after the 
burst. This allowed us to constrain the jet opening angle value. 
We have not noticed differences between the XMM-Newton 
and BeppoSAX samples that can be explained by an absence 
of a selection bias between these observatories. We have not 
observed a significant shrinking of the luminosity distribution 
when we corrected for beaming. This result needs to be con- 
firmed by more burst observations, e.g. by SWIFT. In addition, 
we have found that the SWIFT X-ray afterglows are similar to 
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those of XMM-Newton and BeppoSAX, while Chandra ones 
are biased toward higher luminosities. If this bias is not consid- 
ered, the SWIFT bursts are fainter than the complete sample of 
all bursts listed here. We have also indicated that several dark 
bursts may be simply extinguished in optical due to the high 
absorption observed around the bursts. 
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